Charged current neutrino interactions in core-collapse supernovae in a virial expansion 
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Core-collapse supernovae may depend sensitively on charged current neutrino interactions in 
warm, low density neutron rich matter. A proton in neutron rich matter is more tightly bound 
than is a neutron. This energy shift AU increases the electron energy in v e + n — > p + e, increasing 
the available phase space and absorption cross section. Likewise AU decreases the positron energy 
in i> e + p — > n + e + , decreasing the phase space and cross section. We have calculated AU using a 
model independent virial expansion and we find AU is much larger, at low densities, than the predic- 
tions of many mean field models. Therefore AU could have a significant impact on charged current 
neutrino interactions in supernovae. Preliminary simulations of the accretion phase of core-collapse 
supernovae find that AU increases i> e energies and decreases the v e luminosity. 

PACS numbers: 26.50.+X, 21.65.Ef, 21.65. Cd, 25.30.Pt 



I. INTRODUCTION 

Neutrino and antineutrino capture reactions play a 
crucial role in core collapse supernovae. They provide 
important heating that may reenergize a shock and lead 
to the explosion. In addition, these reactions help de- 
termine the spectra of radiated v e and v e that is impor- 
tant for detection of supernova neutrinos on earth [T] and 
for neutrino oscillations [3] . Finally the rates of neu- 
trino and antineutrino captures determine the ratio of 
neutrons to protons in the neutrino driven wind above 
a proto-neutron star. This is crucial for nucleosynthesis, 
see for example [MB]. 

Supernova neutrinos are radiated from the neutrino- 
sphere, a warm low density gas of neutron rich matter. 
Recently Roberts [7] and Roberts and Reddy [8] have 
suggested that strong interaction energy shifts in this gas 
increase the cross section for 

v e + n^p + e, (1) 

by increasing the available phase space, and decrease the 
cross section for 

D e + p -> n + e + . (2) 

These shifts arise because a proton is more tightly bound 
in neutron rich matter than is a neutron and this binding 
energy difference increases the energy of the e in Eq. ([I]) 
and decrease the energy of the e + in Eq. Martinez- 
Pinedo et al. [9] have also performed supernova simula- 
tions with this energy shift. However, these works [7H5] 
all calculate the shift in a mean field approximation that 
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is very model dependent at low densities. Therefore the 
size of the effect may be poorly determined. 
The energy shift 

AU = U n - U p (3) 

is the difference in potential energy of a neutron U n or 
a proton U p in the medium. This is closely related to 
the symmetry energy that describes how the energy of 
nuclear matter rises as one moves away from equal num- 
bers of neutrons and protons. The symmetry energy has 
been calculated at low densities using a virial expansion 
and found to be large because of correlations and the for- 
mation of bound states such as 4 He, 3 He, and 3 H [lOj . 
Note that these light nuclei may also impact neutrino 
spectra [TT1IT2]. 

Warm low density nuclear matter can be produced in 
the laboratory with heavy ion collisions. Kowalski et 
al. studied collisions of 35 MeV/nucleon 64 Zn projectiles 
with 92 Mo and 197 Au target nuclei and inferred that the 
symmetry energy at densities of 0.01 to 0.05 times satu- 
ration density ranged from 9 to 13.6 MeV |13j . see also 
[14] . These large values are in good agreement with virial 
expansion results. However, they are much larger than 
the predictions of many mean field models, see for exam- 
ple refs. [HHH]. 

In this paper we calculate the energy shifts and in 
medium cross sections for the reactions in Eqs. ([l]),([2]) 
using a virial expansion. The virial expansion makes 
model independent predictions for thermodynamic quan- 
tities based on elastic scattering phase shifts. Some ear- 
lier virial expansion results include the equation of state 
of pure neutron matter |15j , the equation of state of nu- 
clear matter including protons, neutrons, and alpha par- 
ticles [TO], and the long wavelength neutrino response of 
pure neutron matter |16) . Recently, full astrophysical 
equations of state giving the pressure of nuclear matter 
as a function of temperature, density, and proton fraction 
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have been developed that reduce to the virial expansion 
at low densities [T7HI9] . 

This paper is organized as follows. We describe the 
virial expansion formalism in Sec. [TXJ Next, in Sec. |III| 
we present results for the proton fraction in beta equi- 
librium and for the in medium neutrino cross sections. 
Then, in Sec. |IV[ we present preliminary simulations of 
the accretion phase in core collapse supernovae includ- 
ing AU. Finally, in Sec. |Vj we discuss future work and 
conclude. 



In general it is a simple matter to numerically find values 
of z p and z n that reproduce desired n n and n p values in 
Eqs. (|8|,([9|. The chemical potentials then follow from 
Hi = T\m~Zi). Note that \ii do not include the nucleon 
rest mass. 

The virial coefficient b n describes pure neutron matter. 
It is calculated from the observed isospin one nucleon- 
nucleon elastic scattering phase shifts. We fit the nu- 
merical values from ref. |10| over the temperature range 
1 < T < 20 MeV with, 



II. FORMALISM 

We consider a low density warm gas of only neutrons 
and protons. For simplicity we neglect contributions of 
bound states involving three or more nucleons such as 
4 He nuclei. They will be discussed below. Note that 
the effects of deuteron bound states will be implicitly 
included in the virial coefficients, see below. In the virial 
expansion the pressure P is expanded in powers of the 
neutron z n , 



(4) 



and proton z. 



b n (T) 0.3084- 



0.0191 
T 



5.8 x 10~ 4 T. 



(10) 



Here T is in MeV. Likewise the virial coefficient b pn de- 
scribes the interactions between protons and neutrons. 
We again fit the numerical values from ref. [10] . 



b pn (T) 



2 099 

-0.9885 + 2.502 exp[-^=-] - 0.0179T . (11) 



Here the exponential term describes the contributions of 



deuterium bound states. Equations ( 10 ),( 11 ) are good to 
1% over the range 1 MeV < T < 20 MeV. For higher 



temperatures 20 < T < 50 MeV, Eqs. (10), (nil), although 



not explicitly fit to detailed phase shifts, still appear to 
provide reasonable qualitative behavior. 



= e P P /T 



(5) 



fugacities [TU]. Here the neutron chemical potential is //„, 
the proton chemical potential is /i p and T is the temper- 
ature. The virial expansion is valid at low densities and 
or high temperatures where z n , z p < 1 . This roughly cor- 
responds to densities n so that n/T 3 / 2 < 2 x 10" 4 fur 3 
MeV" 3 / 2 . The virial expansion further assumes that a 
phase transition has not taken place from the original 
very high temperature gas phase. 

The pressure, to second order in the fugacities, is [TU] . 

2T 

P = ~j^{ Z n + Z P + ( Z n + z p) b n + 2z p Z n b pn ) , (6) 

where the nucleon thermal wavelength A is 



A 



2tt \l/2 



MT 



(7) 



with M the nucleon mass. We discuss the second or- 
der virial coefficients b n and b pn below. For simplicity we 
neglect the mass difference between the neutron and pro- 
ton. This mass difference will contribute an effect similar 
to (but in general smaller than) the energy shift that we 
consider. The neutron n n and proton n p densities are 

ma. 



A 3 



A 3 



+ 2z 2 b n 



1zpZ n bp n j 



(8) 



(9) 



A. Energy Shift 

To calculate the energy shift we start with the free 
energy density / = {E — TS)/V with E the internal 
energy, S the entropy and V the volume. The use of 
the free energy, instead of the internal energy, will be 
discussed below. In the virial expansion [10] 



/ = n n fi n + n p fJ,p - P . 
Using Eq. ^ we have, 

/ =n n Tlnz n + n p T\nz p 



(12) 



2T 
A 3 " 



{z n + z p + (z 2 + z p )b n + 2z p z n b pn } 



Single particle energies will be calculated from E.- L — 
df/drii. To simplify the calculation we invert Eqs. Q 
and ([9]) to second order in the densities, 



X fir. 



{l - X 3 (n n b n + ripbpn)} 



(13) 



A 3 n 



^{1- \ 3 {n p b n +n n b pn )} , (14) 



giving for the free energy density 



/ w n n Tln 
A 3 T 



n„Tln— : — 
p 2 



- T(n p + n„) 



} {(n 2 n + n 2 p )b n + 2n n n p b pn } +0{n 3 ). (15) 
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This approximation is accurate at very low densities and 
will give us very simple results that provide physical in- 
sight. Later we will obtain more accurate results by ex- 
actly solving Eqs. Q, It is now a simple matter to 
calculate the single particle energies using Eq. (15), 



n X 3 

= Tln-^- X 3 T(n n b n + n p b pn ) , (16) 



Here Gf is the Fermi constant, 9 C the Cabibbo angle, 9 
the scattering angle, and the nucleon axial charge g a = 
1.267. The energy transferred to the medium is q = 
E v — E e for neutrino energy E v and electron energy E e , 
and q is the momentum transferred to the medium q 2 = 
El+E 2 e ~ 2E u E e cos9. Finally, 



f(Ee) = 



1 



Exp[(£ e - /x e )/T] + 1 



(24) 



E„ 



V dn p ) ? 



Tln^ 



X 3 T(n p b n 



(17) 



We measure the energy shift U with respect to the energy 
of a noninteracting Fermi gas. For a free Fermi gas one 
has second virial coefficients ITUl. 



b° =- 



1 

W 2 



b° pn = 0. 



(18) 



is the Fermi Dirac distribution for the outgoing electron. 

Perhaps the simplest model for the response func- 
tion Sv a (qo,q) is to assume the neutrino strikes a heavy 
free nucleon at rest. In this case the response function 
<SW(<70j</) oc 5(qo) so that E e = E v . As we will justify 
below, the effect of AU is to shift the response so that 
Sv m (qo,Q) <H<7o + AU). In this case the energy of the 
outgoing electron will be 



If one expands the pressure of a free Fermi gas in powers 
of the fugacity z, the coefficient of the z 2 term is given 
by b n . Therefore, the single particle energies E® for free 
Fermi gases are given by Eqs. ( 16 ) , ( 171 with b n —> 6° and 
b m 



E P — E v 



AU. 



(25) 



The ratio of the total cross section 
integral of Eq 



angle and energy 
(AU) calculated with AU, to 
the cross section calculated with AU = 0, o v _ (0) will be 



(231, u v 



0. We now have simple results for the neutron and just the ratio of outgoing electron phase spaces, 



proton energy shifts, 



U n = E n - E° = -X 3 T(n n b n + : 



a Ve (AU) (E v + AU) 2 [1 - f(E v + AU)} 



(19) 



,(0) 



^[1-/(^)1 



(26) 



U p = E p -E° p = -X 3 T(n p b n + n n b pn ), (20) 

and finally the difference in energy shifts is 

AU = U n -U p = X 3 T(n n - n p )(b pn - b n ) . (21) 

Here the difference in virial coefficients for interacting 
and free Fermi gases is 



This equation has a simple interpretation. The shift AU 
increases the electron energy and this increases the avail- 
able phase space and therefore the absorption cross sec- 
tion. 

Likewise for antineutrino capture D e +p — > n + e + much 
the same thing happens in reverse. Now the response 
would be approximately (qo, q) cx S(qo — AU) so that 
the positron energy is reduced by the energy shift 



1 

2V2 



E„+ = En - AU . 



(22) 



(27) 



In the next Section we will show that the energy of the 
outgoing electron for v e capture on a neutron will be in- 
creased by AU while the energy of the outgoing positron 
from v e capture on a proton will be decreased by AU. 
Equation (21 1 is a major result of this paper because it 



Therefore the ratio of total cross section cr p<! (AU) with 
AU to the cross section <7 Pe (0) without AU is 



apJAU) = (Eu - AU)' 



0(E P - AU) 



(28) 



provides a model independent prediction for AU in terms 
of virial coefficients calculated from nucleon-nucleon elas- 
tic scattering phase shifts. 



B. Absorption cross sections 

We now calculate the in medium cross section per unit 
volume for v e + n — > p + e [5] . 



cPa 



G 2 cos 2 9, 



V dcos9dE P 



4 _ L , [1 + 3g 2 a + (1 - £)cob0] 

El[l-f(E e )]S v (q ,q) (23) 



Now there is an energy threshold (given by the G func- 
tion) where the cross section is approximately zero until 
En > AU. In Eq. (28) we have neglected Pauli blocking 
for the outgoing positron. Again Eq. (28) has a simple 
interpretation. The shift AU reduces both the energy of 
the outgoing positron and the available phase space and 
this reduces the cross section for antineutrino absorption. 

We now wish to justify the simple results in 
Eqs. (26 1,(28) with a more detailed mean field model of 
the response function S„ B (qo,q). We start with a simple 
model for the in medium single neutron e n (q) and single 
proton e p (q) spectra, 



£n(q) 



2M 



+ U n , 



e P (?) 



T 
2M 



U v 



(29) 
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Here M is the nucleon mass. Note an effective mass M* 



in Eq. (29) is not expected to significantly change our 
results. Furthermore at the low densities that we are in- 
terested in (< O.OOf fin" 3 ), we expect M* w M. There 
is an important consistency requirement between the en- 
ergy shift Ui and the interacting chemical potential /i^. 
The spectrum in Eq. ( 29 1 , with the interacting chemical 



potential, should give the correct nucleon density 



d 3 p 



1 



(2tt)3 Exp(( ei (p) - ih)/T) + 1 



for i = p or n. This requires 



Ui = m 



ri > 



(30) 



(31) 



where fif is the chemical potential of a free Fermi gas. 
Thus the energy shift is just the difference between the 
interacting and free chemical potentials. If Eq. (31) is 



not satisfied, then the mean field response will likely be 
calculated for the wrong nucleon density. Expanding Eq. 
(13) for fx n to second order in the density, and then cal- 
culating by replacing b n — > 6° and b pn — > one has 



Mn = -A 3 T(n n 6 n + n p bp n ) + 0(n 2 ) . 



(32) 



This agrees with Eq. ( 19 ) to lowest order in the density. 



Therefore calculating the energy shift in terms of the free 
energy, see Eq. (15), leads to a consistent definition of 

U^ ^ 

We now calculate the neutrino response S„ c (qo,q) for 
v e + n — > p + e in a mean field approximation assuming 
the spectra in Eq. p9| [20], 



M T 



nq(l 



-In 



e (e mi „-M/T + i 

3(e m i„-An)/T _|_ g-2 



Here, 



is 



M 

v 



(go + AU - 



q 2 n2 
2M> 



}. (33) 



(34) 



and (i n = /ji n — U n . Finally, the detailed balance factor z 
involves the energy transfer and the difference in chemical 
potentials z = (go — [!■„ + fj, p )/T. The response function 
for antincutrino absorption Sp e (qo,q) is obtained from 
Eq. (33) with AU — > —AU and fi n — > fX p , The cross 
section then follows from Eq. (23) with E e replaced by 



the positron energy E e + and f(E e ) replaced by the Fermi 
Dirac distribution for positrons. In the next section we 
will show results for the total cross section obtained by 
integrating Eq. ( 23 1 using Eq. ( [M] ) over outgoing lepton 
energy and scattering angle. These cross sections agree 
closely with the simple phase space ratios in Eqs. (26), 
p). 



III. RESULTS 



In this section we present results for the composition 
of matter in beta equilibrium, neutrino and antineutrino 



absorption cross sections, and for the ratio of cross sec- 
tions calculated with and without energy shifts. We start 
with determining the proton fraction Y p for matter in 
beta equilibrium at baryon density n. The procedure is 
to guess a Y p value and numerically solve Eqs. (13), (14) 



with n p = Y p n and n n = (1 — Y p )n for fi p and \i n . Next 
Y p is adjusted until 



(35) 



Here /j, e is the chemical potential of a relativistic Fermi 
gas of electrons with density n p . 
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FIG. 1: (Color on line) Proton fraction Y p for matter in beta 
equilibrium at a temperature of T — 5 MeV vs baryon density 
n. The black dashed curve is for a free Fermi gas while the 
solid red line includes interactions in a virial expansion. 



In Fig. [T] we show Y p in beta equilibrium versus baryon 
density n for a temperature of T = 5 MeV. This temper- 
ature and range of densities in Fig. [l] are typical for 
the neutrinosphere, at least at early times. We note 
that n = 0.001 fm~ 3 corresponds to a mass density of 
1.7 x 10 12 g/cm 3 . At later times, during protoneutron 
star cooling, the neutrinosphere may move to higher den- 
sities. We discuss this more below. We see that the virial 
interactions increase the equilibrium Y p . 

In Fig. [2] we plot the differential cross section per unit 
volume (l/V)da/dE. This is the angular integral of Eq. 
Q at a density of n = 0.001 fm" 3 , T = 5 MeV, and for 
a neutrino or antineutrino energy of 3T =15 MeV. This 
is a typical neutrino energy for this temperature. Calcu- 



lations that include energy shifts are plotted as thick red 
lines. Note that the energy shifts U are calculated from 
Eq. |3T|). This equation agrees with Eqs. (IT9|), pOl at 



low densities and is more consistent at higher densities 
as we discuss below. 

Cross sections for antineutrinos in Fig. [2] are smaller 
than neutrino cross sections simply because the density 
of protons n p is smaller than the density of neutrons n n , 
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see Fig. [T] Note that the cross section is normalized per 
unit volume rather than per nucleon. This reduction for 
antineutrinos is somewhat mitigated because positrons 
have less Pauli blocking than do electrons. 
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FIG. 2: (Color on line) Differential cross section per unit 
volume for neutrino (solid lines) or antineutrino (dashed lines) 
absorption versus energy E of the outgoing charged lepton. 
The thick red lines include energy shifts and peak at the lowest 
and highest energies. Shifts are neglected for the thin black 
lines that peak near E = 15 MeV. This is for a temperature 
T — 5 MeV, a density of n — 0.001 fm -3 and a neutrino 
(antineutrino) energy of 15 MeV. 

The energy shifts are seen in Fig. [2] to increase the 
energy of the outgoing electron and to increase the cross 
section. Likewise the energy shifts decrease the energy of 
the outgoing positron and reduce the antineutrino cross 
section. We obtain the total cross section as the energy 
integral of da/dE. The ratio of the total cross sections 
with and without energy shifts from Fig. [2] agrees well 
with the simple phase space estimates of Eqs. p6|,(28|. 



The ratio of cross sections with and without energy 
shifts, Eqs. is plotted in Fig. [3] for T = 5 MeV 

and E v = 15 MeV. The proton fraction Y p for all of the 
calculations is the beta equilibrium value in the virial 
expansion (solid red line in Fig. [I]). The thin lines use 
the lowest order energy shifts from Eqs. (19 1,(20) while 



the thick lines use the difference in chemical potentials 
from Eq. (31 1. These two prescriptions for the energy 
shifts agree at very low densities. However, at higher 
densities the lowest order energy shifts predict a neutrino 
ratio that is larger than that for Eq. (31 1. Since this later 



energy shift insures that the density is correct, we believe 
the thick lines in Fig.|3]are the correct predictions for the 
neutrino and antineutrino ratios. 

At a density of n = 0.001 fm~ 3 (1.7 x 10 12 g/cm 3 ) the 
neutrino absorption cross section is enhanced by 73% be- 
cause of the energy shifts, while the antineutrino absorp- 
tion cross section is reduced by 28%. We now compare 




v + p — > n + e 
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FIG. 3: (Color on line) Ratio of absorption cross sections with 
and without energy shift versus density n. The solid black 
lines are for neutrinos, Eq. (26 1, while the dashed red lines 
are for antineutrinos, Eq. (281. The thick lines use energy 
shifts from Eq. (311 while the thin lines use en ergy s hift s 
valid to lowest order in the density from Eqs. (19 1, (20 1. 



Finally, the dot-dashed (neutrino) and dotted (antineutrino) 
curves use the small mean field energy shifts from Eq. ( |36[ |. 
This is for T = 5 MeV and for a neutrino energy of 15 MeV. 



these virial expansion results to a mean field model. For 
example, Roberts and Reddy [5] consider a mean field 
model GM3 that is normalized to nuclear phenomenol- 
ogy at saturation density n Q — 0.16 fm~ 3 . In this model 
the energy shift AU m f is 



AU mf = 40- 



MeV. 



"o 



(36) 



This model, compared to the virial expansion, predicts 
much smaller energy shifts. Furthermore, the ratios of 
neutrino or antineutrino cross sections in Fig. [3] are much 
closer to one . For example at n — 0.001 fm , neutrino 
absorption is only enhanced by 6%. Alternatively, the 
mean field model IUFSU [53] has a nonlinear density de- 
pendence and predicts an energy shift larger than that for 
GM3 but still less than the virial predictions. In Table [I] 
we collect energy shifts for n — 0.001 fm -3 and T = 5 
MeV. 



Model 


AU (MeV) 


Lowest order virial, Eq. (21 1 


3.85 


Virial lH-l4, Eq. (311 


2.27 


Mean field model GM3, Eq. (36 » 


0.23 


Mean field model IUFSU [24] 


1.11 



TABLE I: Energy shift AU predicted by different approaches 
at a density n = 0.001 fm -3 and a temperature T — 5 MeV. 
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The virial expansion calculates the pressure as a power 
series in the fugacities that is valid at low densities and 
or high temperatures. For the conditions in Fig. [3j the 
neutron fugacity is z n < 0.16 while the proton fugacity 
is z p < 0.0045. These small values suggest that higher 
order corrections to the virial expansion will be small. 

The inclusion of alpha particles was shown in ref. [TU] 
to significantly enhance the symmetry free energy at 
higher densities and or lower temperatures. Therefore we 
expect the formation of alpha particles (and other nuclei) 
to further enhance the virial energy shifts and neutrino 
absorption cross section changes. Although preliminary 
calculations find small effects for the conditions in Fig. [3j 
alpha particles do make significant contributions at lower 
temperatures and or higher densities. 

In this paper we have only calculated the effects of 
an energy shift AU on neutrino interactions. In addi- 
tion, strong interactions will introduce other correlations 
between nucleons and these will impact neutrino interac- 
tions. Often these correlations are included in model de- 
pendent RPA calculations, see for example [2"TM23| . How- 
ever the effects of correlations have been calculated for 
the long wavelength neutrino response of pure neutron 
matter using the virial expansion [16] . In future work we 
will calculate the effects of correlations in nuclear matter 
using the virial expansion for both charged current and 
neutral current interactions. 



IV. SUPERNOVA SIMULATIONS 

We perform exploratory simulations to gauge the influ- 
ence of the energy shift on the neutrino signal during the 
accretion phase of core-collapse supernovae. We make 
use of nuGRlD fZT\ [55], a spherically-symmetric, general- 
relativistic, Eulerian hydrodynamics code with a two- 
moment neutrino radiation transport solver. For these 
simulations we take the standard 15 Mq, solar metallic- 
ity, core-collapse supernova progenitor profile from |29j . 
We employ the Lattimer and Swesty [3D] equation of 
state with an incompressibility modulus of 220 MeV. The 
neutrino interaction rates are generated using NuLib, 
an open-source neutrino interaction library available at 
http://www.nulib.org. In NuLib, the absorption cross 
sections for neutrino and antineutrino capture on free 
neutron and protons are taken from ref. [31J . 

<l - °^(l + 39l)El (l~|^) 1/2 

xW M [l- f(E e -)}, (37) 

and 

< = V (1 + 3 ^ + ( 1 -||) 1/2 

xW M [l-f(E e+ )], (38) 



where Wm is the weak magnetism correction [5J, and 
f(E e -) and f(E e +) are the Fermi Dirac distributions for 
the outgoing electron and positron, respectively. We set 
E e - =E V + A np + AU and E e + = E v - A np - AU. Here, 
A np is the neutron-proton mass difference and AU is the 
energy shift. For simplicity in the implementation, and 
to remain conservative in the calculation of the energy 
shift, we use the following definition for AU, which fol- 
lows closely from Eq. ( ]21[ ) . 

AU = T Min(A(T) 3 (n„ - n p ), 1) [b pn (T*) - b n (T*)} , 

(39) 

where the Min function is to limit the energy shift in re- 
gions where the virial approximation breaks down. Like- 
wise, we use T* = Max(T, 1 MeV) to limit the size of the 
virial coefficients for low temperatures. The correspond- 
ing neutrino emissivity from electron/positron capture 
on free nucleons is calculated consistently using detailed 
balance [281I3T]. 
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FIG. 4: (Color on-line) Average energy (left panel) and lumi- 
nosity (right panel) of electron-type neutrinos and antineutri- 
nos with and without energy shift. The solid black lines are 
for neutrinos, while the dashed red lines are for antineutrinos. 
The thick lines are from the simulation including the energy 
shift, while the thin lines do not include the energy shift. 

In Fig. |4j we show the influence of the energy shift 
on the preexplosion neutrino signal through the evolu- 
tion of the average neutrino energy (left) and neutrino 
luminosity (right) for both the electron neutrinos (solid, 
black lines) and antineutrinos (dashed, red lines). The 
thick (thin) lines are the results with (without) the en- 
ergy shift. The third species included in the simula- 
tions, a characteristic heavy-lepton neutrino, shows very 
little change with the inclusion of the energy shift and 
is not shown. The energy shift leads to an increase in 
the average energy of the electron antineutrinos. This 
is a result of a decrease in the absorption opacity of the 
electron antineutrons on free protons, moving the neutri- 
nosphere for a given neutrino energy to lower radii and 
therefore higher matter temperatures. The magnitude of 
this difference is ~0.25MeV, (~ 1.8%) at 100 ms after 
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bounce. The difference grows as the neutrinosphere re- 
cedes to higher densities where the energy shift is large. 
After 450 ms of postbounce evolution, the difference is 
^0.94 Me V (~ 5%). Note that neutrino-electron scatter- 
ing can lower the v e average energy and could somewhat 
limit the impact of AU on the D e spectrum. The av- 
erage energy of electron neutrinos shows no dependence 
on the energy shift. The high free neutron fraction in 
the postshock region forces the electron neutrino neutri- 
nospheres to lower densities where the influence of the 
energy shift is not as effective. The other strong effect 
of the energy shift is the reduction of the electron neu- 
trino luminosity. The decrease is 1.5 Bs -1 (~ 3%) at 
100 ms after bounce and increases to 2.5 Bs -1 (~ 10%) 
at 450 ms after bounce. We attribute the reduction to the 
lower electron fraction found throughout the postshock 
region, but most importantly around the neutrinospheres 
(where AY e /Y e between the two simulations is ~ 20%). 
The lower electron fraction is a result of the overall lower 
electron antineutrino opacities and emissivities when the 
energy shift is included. We note that most (~90%) of 
the electron-neutrino luminosity change already obtains 
when including the shift only in the antineutrino opaci- 
ties. 

The potentially significant quantitative changes that 
arise as a result of the energy shift, both during the pre- 
explosion phase (as shown here) and presumably during 
the postexplosion phase warrant a much more in-depth 
analysis in future work. 



V. DISCUSSION AND CONCLUSIONS 

In this section we discuss our results and conclude. 
First, as shown in Fig. [3] the virial expansion energy shift 
AU, at low densities, is much larger than that predicted 
by many mean field models and leads to much larger 
changes in neutrino cross sections. This is a major con- 
clusion of the present paper. Furthermore the virial pre- 
dictions are based on elastic scattering phase shift data 
and are model independent. Our preliminary simulations 
of the accretion phase of core collapse supernovae find 
that AU increases the D e energies and decreases the v e 
luminosity, see Fig. [4] 

We expect the energy shifts to have even larger effects 
on neutrino interactions at higher densities n > 0.001 
fm~ 3 . However the virial expansion itself may not be 
directly applicable at these densities. Therefore, it is im- 
portant to calculate the properties of warm neutron rich 
matter in other microscopic approaches. Although con- 
ventional quantum Monte Carlo approaches, such as ref. 
[25|, I26j , may have difficulties with both the nonuniform 
matter and the high temperatures, other techniques may 
be more promising. For example lattice effective field the- 
ory (32J, [33] should be directly applicable. Furthermore, 



tentials in a microscopic approach in order to determine 
the energy shifts and their impact on neutrino absorp- 
tion. 

We emphasize that matter at neutrinosphere temper- 
atures and sub-saturation densities can be directly pro- 
duced in the laboratory with heavy ion collisions, see for 
example ref. |13j . Furthermore, new radioactive beam 
facilities will allow the study of more neutron rich con- 
ditions. Terrestrial experiments can probe the equation 
of state, symmetry energy, and composition of neutrino 
sphere like matter. These properties arc important for 
neutrino interactions in core-collapse supernovae. 

Large energy shift effects could lower the energy or 
luminosity of v e radiated in supernovae. Although about 
20 v e were detected from SN1987a, we have almost no 
experimental information on v e energies from SN1987a. 
Therefore it is important to have a supernova v e detector 
with good energy resolution (such as liquid Ar 134,135]) to 
complement Super Kamiokande and other existing good 
v e detectors. We also note that Pb based detectors such 
as HALO [36] have v e sensitivity, while the energy of 
Ufj, and v T can be measured with neutrino-nucleus elastic 
scattering detectors [37] for example. 

A reduction in the v e energies or luminosity will likely 
make the neutrino driven wind above a protoneutron star 
more neutron rich. This is important for nucleosynthesis. 
Perhaps energy shift effects are large enough to provide 
the necessary free neutrons in order for the r-process to 
occur in the neutrino driven wind. Our preliminary simu- 
lations of the accretion phase of core collapse supernovae 
should be extended to the explosion and protoneutron 
star cooling phases using energy shifts that are calcu- 
lated accurately at higher densities where the virial ex- 
pansion is not directly valid. We will present additional 
simulation results in a later publication. 

In conclusion, a proton in neutron rich matter is more 
tightly bound than is a neutron. This energy shift AU in- 
creases the electron energy in v e + n — > p + e, increasing 
the available phase space and absorption cross section, 
see Eq. (26 1. Likewise AU decreases the positron en- 
ergy in v e + p — > n + e + , decreasing the phase space and 



cross section, see Eq. (28). We have calculated AU using 



Eq. (31 ) shows that the energy shifts can be directly de- 
termined from chemical potentials. Therefore one only 
needs to calculate the neutron and proton chemical po- 



a model independent virial expansion and we find AU 
is much larger, at low densities, than the predictions of 
many mean field models. Therefore AU could have a sig- 
nificant impact on charged current neutrino interactions 
in supernovae. 
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